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1. Introduction 

Spacetimes admitting naked singularities provide test-beds for the cosmic censorship 
conjecture. Of course the mere existence of an example of a spacetime admitting a 
naked singularity (NS) is not evidence that the conjecture is invalid, as such spacetimes 
are typically unrealistic in that they possess a high degree of symmetry. Rather, the 
conjecture may be probed by determining whether or not these spacetimes (and their 
NS) are stable under perturbations. This addresses the question of whether or not NS 
may arise from open sets of initial data, which in turn addresses (to an extent) the 
question of whether or not NS may arise in nature. For example, the NS of a charged 
spherical black hole cannot arise in nature as it has been shown that the Cauchy horizon 
accompanying the singularity is unstable under small perturbations (see jT] for a review). 
However, there is evidence that this instability is not present for certain spherically 
symmetric self-similar spacetimes. 

Harada and Maeda have shown the existence of perfect fluid spacetimes with a soft 
equation of state that admit naked singularities but which are stable under perturbations 
impinging on the singularity: these have been dubbed the general relativistic Larson- 
Penston spacetimes (GRLP) [2]. (The fate of perturbations impinging on the Cauchy 
horizon has not been studied, and this is an important question.) It has also been 
shown that individual massless scalar wave modes of the form e''^'"(j){x) remain finite in 
the approach to the Cauchy horizon OC ^ Xf^ 111 cL wide class of spherically symmetric 
self-similar spacetimes Here, v is an advanced Bondi coordinate and x = v/r is 
a similarity variable (see below; r is the radius function of the spherically symmetric 
spacetime). The same has also been shown to be true for the individual modes of 
gravitational (metric and matter) perturbations of the self-similar Vaidya spacetime |3]. 

The aim of the present paper is to expand upon the second of these three points 
by studying the propagation of a minimally coupled, massless scalar field in self-similar 
spacetimes without recourse to a mode decomposition (Fourier transform). We find that 
the results for individual modes essentially carry over to the multipoles of the full field, 
and so have greater significance. (We us the term multipole to refer to the coefficients 
(pi^m of the spherical harmonics in the standard angular decomposition of a scalar field 
on a spherically symmetric background: (p = Ylem^i,m^e^i^^v)-) "^^^ norm of 
each multipole field, its pointwise values and its local energy density remain finite in 
the approach to the Cauchy horizon. We consider this to be evidence of the stability of 
such naked singularities: the scalar field can be considered to model the behaviour of 
gravitational perturbations of the spacetime. 

In the following section, we describe the global structure of the class of spacetimes 
under consideration: spherically symmetric spacetimes admitting a naked singularity 
and whose energy momentum tensor obeys the weak energy condition. In section 3, 
we discuss the and local stability of a massless scalar field propagating in such a 
spacetime, and study the local energy density of the field in section 4. We consider 
massive fields in section 5, and discuss our results in section 6. We use the curvature 
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conventions of [S], setting G = c = 1, and use the notation of ||6| for Sobolev spaces. A 
black square ■ indicates the end or absence of a proof. 

2. Self-similar spherically symmetric spacetimes admitting a naked 
singularity. 

We will consider the class of spacetimes which have the following properties. Spacetime 
{Ai,g) is spherically symmetric and admits a homothetic Killing vector field. These 
symmetries pick out a scaling origin O on the central world-line r = (which we 
will refer to as the axis), where r is the radius function of the spacetime. We assume 
regularity of the axis to the past of O and of the past null cone Af of O. We will use 
advanced Bondi co-ordinates [v, r) where v labels the past null cones of r = and is 
taken to increase into the future. Translation freedom in v allows us to situate the 
scaling origin at (f,r) = (0,0) and identifies v = with J\f. The homothetic Killing 
field is 



where dQ"^ is the line element of the unit 2-sphere. The homothetic symmetry implies 
that G{v,r) = G{x),ip{v,r) = ip{x) where x = v/r. The only co-ordinate freedom 
remaining in is then v V{v); this is removed by taking v to measure proper time 
along the regular center r = 0. 

We will not specify the energy-momentum tensor of {A4,g), but will demand that 
it satisfies the null energy condition: Tabk°'k^ > for all null vectors k. Note that this 
condition is implied by both the strong energy condition and the weak energy condition. 
A complete description of energy conditions in spherical symmetry is given in [3]. Of 
these, the following will be used. 



We define the interior region Aimt of spacetime to be the interior of A/", i.e. the 
interior of the causal past of O. We demand that there are no trapped 2-spheres in 
A4int- This is an initial regularity condition, and is equivalent to G > in Aimt- We 
also note that except for the trivial case of fiat spacetime, there must be a curvature 
singularity at C jH]. 



The presence of a naked singularity is characterised by the following result proven 
in[al. 




(1) 



xip' < 0, 

G' -^'{l-xG)G < 0. 



(2) 

(3) 



Theorem 1 There exists a future pointing radial null geodesic in Ai with past endpoint 
on O if and only if there exists a positive root of the equation xG{x) — 1 = 0. ■ 
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The features above define the class of spacetimes under consideration in the 
remainder of this paper: {Ai,g) is spherically symmetric and self-similar with line 
element (^. There exists a naked singularity, so that there is a smallest positive 
root Xc of xG — 1 = 0. Then the Cauchy horizon Ti. is x = Xc- We take G,i/j 
to be analytic on (—00, Xc]. Einstein's equation and the null energy condition are 
assumed, so that Q and Q apply. The corresponding conformal diagram is shown 
in Figure One. There are subsets of the following categories of spacetimes included in 
this class: (i) Vaidya spacetimes (spherically symmetric null dust); (ii) Einstein-Klein- 
Gordon spacetimes (spherically symmetric massless scalar field coupled to gravity); 
(iii) Lemaitre-Tolman-Bondi spacetimes (spherically symmetric dust); (iv) perfect fluid 
filled spacetimes (where self-similarity enforces a linear equation of state - see [7j); 
(v) Einstein-S'[/(2) spacetimes (spherically symmetric SU{2) sigma model coupled to 
gravity - see |H]). Of these, the metric functions G^ip are explicitly available only for 
Vaidya spacetime, which has = 0, 2G = 1 — 2Ax, where A is a constant. A naked 
singularity occurs for A G (0,1/16). Among the other cases, the similarity coordinate 
z = t/r (where t is proper time along the fluid flow lines in the dust and perfect fluid 
cases and a polar slicing time coordinate in the sigma model case) is more natural 
than our x, and a nontrivial coordinate transformation is required to determine if), G 
explicitly. We emphasise however that this explicit representation is not required for 
the present purpose. 



r = 0, u > 




Figure 1. Conformal diagram for an example of a self-similar spacetime admitting a 
globally naked singularity. We use the advanced Bondi co-ordinates v and r described 
in Section 2. The Cauchy horizon is shown dashed, the event horizon as a double line 
and the apparent horizon as a bold curve. M is the past null cone of the scaling origin. 
Other structures can arise; there may be no apparent or event horizon; the censored 
portion of the singularity may be null; the naked portion of the singularity may be 
time- like. 



The temporal nature of the surfaces x = constant is central to the discussion below, 
and so we note the following which is immediate from (^. 
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Lemma 1 The surface x = Xq is space-like (respectively time-like, null) if and only if 
Xq{1 — XqG{xq)) > (respectively, < 0, = 0). ■ 

The definition of Xc tlien yields tlie following: 
Corollary 1 The surfaces x =constant are space-like for x G (0, Xc). ■ 

It is clear that at the past null cone x = we have xG — 1 < 0, and so a naked 
singularity forms when the function G meets 1/x from below. We can then show that 
1 + x^G'^Xc) > 0. The case of equality appears to be quite special: by appealing to the 
field equations, it can be ruled out when the energy momentum tensor is that of dust 
(Lemaitre-Tolman-Bondi spacetime), null dust (Vaidya spacetime), a massless scalar 
field or an SU{2) sigma model. (The case of perfect fluid is more complicated, but by 
appealing to a degrees of freedom argument, it appears that equality is non-generic.) 
We will therefore assume that equality does not hold - i.e. we will impose the condition 

1 + xlG'ixc) > (4) 

as one of the defining characteristics of the class of spacetimes being studied. 

Finally, we give a result from [Sj that plays an important role in the derivation of 
the stability results below. 

Lemma 2 G" < prior to the formation of a Cauchy horizon and hence G'{xc) < 0. ■ 

We note that if G'{xc) = 0, then Rabk"'k''\^ = 0, where k'^ is tangent to the 
outgoing radial null direction. This implies that there is no ingoing radiative flux of 
energy-momentum crossing the Cauchy horizon. We rule out this situation as being 
physically unrealistic and so we will assume that 

G'{x,) < 0. (5) 

Note also that this implies that if G{xi) = for some xi > 0, then G{x) < for all 
X > Xi and a naked singularity cannot arise. Thus we must have G{x) > in [0,Xc]. 

3. and local bounds for the scalar field 

Before proceeding, we make some comments on the analysis that is carried out below. 
This analysis is based on the use of naturally arising space-like slices of the form x = 
constant, where x is the homothetic coordinate mentioned above. These surfaces meet 
the naked singularity at its past endpoint {v,r) = (0,0) (the scaling origin). We study 
the evolution in x of the L^— norm of the multipoles of a minimally coupled scalar 
field 0, where the norm is taken with respect to the measure dp = r~^dr. (We deal 
principally with the massless case, and comment later on the massive case.) Note then 
that 4>\x=xo ^ i^^(ffi) requires vanishing of in the limit as r — on the slice x = Xq. 
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This is an undesirable feature, as this constrains to have support outside the past null 
cone of the origin. Wc would much prefer to be able to deal with the case where is 
non-zero on and inside the past null cone. To get around this problem, we introduce 
the field = (pe'^'', k, > and show that maintains finite norm in the evolution, at 
least for sufficiently small values of k. This allows (p to be non-zero at r = 0. 

Our strategy runs as follows. By rewriting the wave equation governing the 
evolution of the field in first order symmetric hyperbolic form, we can apply a standard 
theorem to obtain existence and uniqueness on [xi, Xc) x R where x — Xi & (0, Xc) is an 
initial slice. Applying another standard result, we obtain an a priori bound for a certain 
energy norm, but with a coefficient that grows exponentially as the Cauchy horizon Ti 
is approached. However, in a neighbourhood of 7i, we can introduce a different energy 
norm whose growth is bounded. From this, we prove that has finite H^'^ norm at the 
Cauchy horizon, i.e. 



/ \4>{^c,p)\'^ + \4>A^c,p)\'^ dp <oo. 



A Sobolcv-type inequality then immediately proves boundcdncss of (f){xc, p) for all 
p = Inr G M. These results apply for all solutions generated by smooth data with 
compact support. By taking a sequence of such solutions, we can generalise to the case 
where the initial data - and corresponding solutions - have 0la;=constant ^ -ff^'^(M). 
This allows that is non-zero on the past null cone. 

In the remainder of this paper, we use as coordinates x — v/r and p = Inr. The 
line element then reads 

ds^ = e'^f[-2Ge'l'dx'^ + 2e^{l - 2xG)dxdp + 2xe^(l - xG)dp^ + d^f]. (6) 

We consider the evolution of a massless, minimally coupled scalar field on the spacetime 
with this fine element, restricting attention to the region M. — {A4 : < x < Xc}. Note 
that by Corollary 1, x is a time coordinate in this region. We decompose the scalar field 
into angular modes: = Re^^^ 0^^(x, p)Y^^(^, </?). The angular mode indices i,m 
will be suppressed henceforth, and we take to be the real part of 0^^. Then satisfies 
the 2-dimensional hyperbolic equation 

a0,,, + 2/30,,^ + 70,^^ + {a' + 2/3)0,, + (/3' + 27)0,^ - i{i + l)e^0 = 0(7) 

and = (pe'^f satisfies 

a0,,, + 2/30,,p + 70,pp + {a' + 2(1 - k)/3)0,, + (/?' + 2(1 - 

+ [{k^ - 2K)'y - K(3' - £{£ + l)e^] (P = (8) 



where 



a{x) ^ -2x{l- xG), (9) 
(3 = (3(x) = 1 - 2xG, (10) 
7 = 7(x) = 2G. (11) 
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Here and throughout, a comma denotes a partial derivative. The fact that the 
coefficients are independent of p is a consequence of self-similarity. We also note the 
important fact that a < for x G (0,Xc). 

The most notable feature of (jHl) - which includes ((Tj) as the special case k = - 
is that it is singular at Ti: the space-like surfaces x = constant become characteristic 
in the limit x — > Xc- This is the fundamental difficulty in dealing with the equation. 
In order to obtain an analytic system, we rescale the time coordinate x, placing the 
Cauchy horizon at infinity. This introduces the possibility that the energy will have an 
infinite amount of time to grow, and so will diverge at the Cauchy horizon. However, 
this turns out not to be the case. 

Let Xi G (0, Xc) and define 

(12) 



t 



a[s\ 



By the definitions above, t is an analytic function of x on [xj, Xc), t{xi) = and t oo 



as X 



Let 



\ 
I 



(13) 



Then (jH} may be written in the ffist order symmetric hyperbolic form 

+ A^^p + B0=O, (14) 

where A,B are smooth bounded matrix functions of t on (0, oo), and A is symmetric 
with real, distinct eigenvalues: 

/ \ 

-(3 1 



A 



B 



1 - 







V 



-a{{n^ - 2k)7 - K(3' - i{i + l)e^) 




1 -P 
-2{1-k)/3 2(1 -k) 




The following theorem underpins our subsequent work, and is a standard theorem for 
symmetric hyperbolic linear systems (see e.g. Chapter 12 of |9,). We recall that for Q 
a connected subset of MJ^, C^{Q,X) is the set of smooth functions from Q to X which 
vanish outside some compact subset of Q. (When X = M, it will be omitted.) 



Theorem 2 Let f,g E C^(]R, M'^). Then there exists a unique solution E C°°(M x 
[0,oo),]R^) of the initial value problem consisting of (|i^| ) with the initial conditions 
ip\t=o = f,0,t\t=o = 9- For each t > 0, the function 'fit,-) : M ^ has compact 
support. ■ 
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Corollary 2 Let Xi G (0,Xc) and let f,g & C^(M, M). Then there exists a unique 
solution (f) G C°°(M x [xj, Xc), K) of the initial value problem consisting of (0j with the 
initial conditions 4'\x=Xi = f,4',x\x=Xi = 9- For each fixed x G [xi,Xc), the function 
■) : M — s> M has compact support. ■ 



This corollary is worth stating as it provides, in terms of the natural coordinates, 
the basic existence and uniqueness theorem for the field (f) on the region bounded to the 
past by the past null cone and to the future by the Cauchy horizon. 

The following corollary provides initial control over the growth of the norm of 
0, but note that it provides no information (other than subexponential growth) about 
the norm in the limit as the Cauchy horizon is approached. We define 

m = E[0]{t) = [ M'dp. (15) 

The proof of this corollary is standard, relying on the symmetry of A and the 
boundedness of B in ()14|1 and the fact that the solution of Theorem 1 has compact 
support on each t =constant slice. See for example Theorem 1 of Chapter 12 of jH]. In 
this corollary and throughout, the symbols Co, Ci, . . . will be used to represent (possibly 
different) positive constants that depend only on the metric functions ip, G and on the 
angular mode number i. 

Corollary 3 Subject to the hypotheses of Theorem{^ E[0]{t) is defined for all t > 
and satisfies 

^[^](t)<e^«*E[^](0), 
where Bq = sup^^g I " 25 1 < +oo. Consequently, 

[ \4>{x,p)\'dp <e^°*E(0), 
Jr 

[ |0>,p)pdp<e^"*E(O), 
Jr 

[ |0,,(x,p)|2cip<Cie^°*E(O). 
Jr 



We note that Bq is finite by virtue of the fact that (3, tjj are analytic and bounded 
on (0, Xc] and that i, k, are finite. The bounds on and its derivatives come straight 
from the definition of E{t): the third requires the use of Minkowski's inequality and 
incorporates the exponential growth of as t oo. 

Our next step is to introduce a different energy norm, which can be shown to display 
only bounded growth in the approach to the Cauchy horizon. For this, we return to the 
original second order form (jH)) of the equation for (p. For the solutions of Corollary |21 
the following integral is defined and is differentiable for all x G [xi,Xc)'- 

E{x) = E[4>\{x) = [ -a4>l + -f4>l + dp, (16) 
Jr 
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where H = i{i + l)e'^ — (k^ — 2/t)7 + From the definitions of f3 and 7, we have 

H = i{i+ l)e^ + 2k{1 - k)G - 2/txG', 

and so if < K < 1, then E{x) is non-negative for all x G [xj, Xc). The growth of E{x) 
in the approach to the Cauchy horizon is controlled by the following result. 

Lemma 3 Let 



and let 



< K < min{l, — ^} 



I -2k -2x,{kG'{xc)+x,G"{xc)). 

mo > maxjO, , }. 

Xc 1 - X^G'(Xc) 



Then there exists Xi G (xj,Xc) sttc/i that i?[0](x) zs non-negative on [xj,Xc) and satisfies 

^ < mo^ 
ax 

/or a// X G (xi, Xc). 

Proof: Non-negativity of follows immediately from its definition and from k G [0, 1). 
Smoothness of means that E is differentiable, and that its derivative may be calculated 
by differentiation under the integral. This is simplified in three steps: (i) Integration 
by parts of the term 0,p0,pa; and the removal of a boundary term - permitted as has 
compact support on each slice x = constant, (ii) Removal of the term with (p^^x by 
application of the equation (jHJ. (iii) Removal of a total derivative containing 4'^x4',px- 
This results in ^ = J^I dp, where 

I={a' + 4(1 - k)(3)^% + 2{(3' + 2(1 - k)7)0,x0,p + 
+ [£(^ + l)^'e^ - 2k{kG' + xG")]0^. 

For any tuq G M, we define //? by / = tuqIe + Ir, where 



p 



Ie = -a4>% + ^4>l + {£{£ + l)e^ - {k^ - 2/€)7 + k^Y^ 



so that E = J^lEdp. The coefficient of 0^ in Ir is T(x) = £(£ +1) (■?/'' — mo) e'^ — 2/tTi(x) 
where 

Ti = (1 — K)mQG + (k — mQx)G' + xG". 
The assumed bounds on k and mo yield the following: 

Ti{xc) = (1 - k)— + (k - moXc)G''(xc) + XcG"{xc) 

Xr 



'l-xlG'{Xr 

> 



2Xr 



mo + kG\xc) + XcG"{xc) > 0. 



Hence by continuity, there exists X2 G [xi,Xc) such that Ti(x) > for all x G (x2,Xc]. 
Then using the energy condition ^ and mo > 0, we find T(x) < on [x2, Xc] and so on 
this interval, 

Ir < («' + moa + 4(1 - k)P)^1 + 2{(3' + 2(1 - k)^)^^,^ + (7' - mo7)0'p. 
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Consider the 1-parameter (x) family of quadratic forms in R, S defined by 

Q{R, S; x) = {a' + moa + 4(1 - K)i3)R'^ 

+ 2(/5' + 2(1 - k)^)RS + (7' - mo7)5'. 

The assumption of differentiabihty of G{x) at the Cauchy horizon means that the 
coefficients of Q{R,S]x) are continuous functions of a; at a: = Xc, and using (p iTT|) 
we can calculate 

Q{R, S; Xc) = 2{2k - 1 + a;^G"(x,))i?^ 

+ —(1 -2k- xlG'{xc))RS + 2{G'{xc) - —)S\ 

Xq Xc 

The discriminant of this quadratic form is given by 
16 

A = —(1 -2k- xlG'{xc)){l -2k- rriQXc). 

Now let < K < ^'^''f ^^'-^ (recall that G"(x,) < 0). If 5 = 0, then 

Q{R, 0; Xc) = 2{2k - 1 + xlG\xc))R^ < 0, 

with equality only if i? = 0. If S* 7^ 0, then the calculation of the discriminant 
above shows that Q{R, S; Xc) is negative definite, provided 1 — 2k — rriQXc < 0, i.e. 
if mo > (1 — 2k)/xc- With these inequalities in place, we see that Q(-R, S; Xc) is negative 
definite. Hence by continuity, there exists X3 G [xi^Xc) such that Q{R,S;x) is negative 
definite for all x G [xi,Xc). Hence Ir < on [xi,Xc), where xi = max{a;2,a;3} giving 
/ < rriolE, and the lemma is proven. ■ 

We can now prove the first of the main results of this section. 

Theorem 3 Let cf) he a solution of (0j that is subject to the hypotheses of Corollary{^ 
and of Lemma\^ Then the energy E[(j)]{x) of (p satisfies the a priori bound 

E[4)]{x) <GiE[0]{O), xe[x,,Xc\. 

In particular, has finite H^'^ norm at the Cauchy horizon: 

I ^\xc, p) + ~4>l{xc, p) dp < G2E['?m < 00. (17) 

Proof: From Corollary 13 we obtain 

E[^]{x) < c(a;)E(0), 

where c is a smooth positive function on [xj, Xc) that diverges in the limit x Xc- Then 
< c(xi) < 00, where xi is the value of x identified in Lemma El For x G [xi.Xc), we 
can integrate the result of this lemma to obtain 

E{x) < e^'^^^^Eixi). 

Combining these two inequalities yields 

E[4>]{x) <GE[0]{O), xe[xi,Xc), 
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where C = e'^'^^^''~^'^^c{xi). Since the right hand side of the last inequahty is independent 
of X, we can take the hmit x Xc on the left to see that it applies for all x G [xj, xd- 
The inequality ()17|1 then follows from the definition of E[(f)]{x). ■ 

To obtain a pointwise bound for ^It^, we recall the following Sobolev inequality (see 
p. 1057 of [ID]). 

Lemma 4 Let u E Co~(M). Then for all s eR, 

\u{s)\^<-{ [ \u{t)\Ut+ [ \u'{t)\^dt\. (18) 
2 Uk Jr J 



Theorem 4 Let (j) be a solution of (0j that is subject to the hypotheses of Corollary\^ 
and of Lemma{^ Then (f) is bounded on [xjjXc] x R. 

Proof: Theorem ini shows that has finite H^'^ norm on each slice [xi,Xc) 3 x = 
constant. Indeed we can write 

Jr 

for some constant C3 > that depends only on G, if) and Applying Lemma ID yields 
10(2^5 p)| < ^-^(0), X E [xj,Xc),p E M. As the bounding term is independent of x, we 
can take the limit x ^ Xc and so extend the bound to the Cauchy horizon x = Xc- ■ 

To conclude this section, we point out that the results of Theorem El and Theorem 
m can be extended to solutions that lie in if^'^(M) on each slice x = constant. For 
these solutions, the field need not vanish at the origin. 

Theorem 5 Let f E H^'^{R) and g E ^^(M). Then for all k E [0,min{2, (1 - 
x^G"(xc))/2)}] the initial value problem consisting of (0) with data (t)\x=x^ = f, (p,x\x=Xi = 
g has a unique solution (p E C([xj, xj, -ff^'^(M)). For such solutions, the field values 
(j){x,p) and energy E[(j)]{x) satisfy a priori bounds 

4>{x,p) < CiE[0]{O), {x,p) E [ 
E[4>]{x) <C2E[0]iO), xE[x,,x,]. 

Proof: The space C^(M) is dense in both H^''^(R) and L^(M) with the appropriate 
norms (see Corollary 2.30 of [U] for the former and Theorem 3.17 and Corollary 3.23 of 
jH] for the latter result). Thus we can take a Cauchy sequence {/(m)}m=05 (respectively 
{g(m)}m=o) of functions in C^(M) which converges in the H^'^ norm to / (respectively, 
in the L^(M) norm to g). For each m > 0, the functions f{m),g{m) can be taken as initial 
data for 0, respectively for the equation (jH}. Then the hypotheses of Corollary 2 and 
Theorems 3 and 4 above are satisfied, and we obtain a sequence of solutions {0(m)}m=o 
with 0(m) E C°°{[xi,Xc) X M, M) for all m > 0, and such that each 0(m) has compact 
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support on each x = constant slice. Furthermore, the bounds of Theorem El apply for 
each m: 

E[4>(m)]{x) < Ci^[(^(„)](0), X G [Xi,Xc]. (19) 

Note that the constant Ci is the same for each m, and so by linearity 

)](0), X G [XjjXj, 

for all m, n > 0. The convergence properties of {/(m)}m=0' ^^"^ {5'(m)}m=o imply that 
the sequence of real numbers {E[(p(^rn)\iS^)}m=Q converges, and so for each x G [xi,x^, 
0(m),x)}m=o a Cauchy sequence in the norm 

II (m, v)\\e= [ -av^ + lul + + l)e'^ - {k^ - 2/^)7 + k(5')u\ 
Jm. 

defined so that \\{^,^,x)\\e = E[(f)]{x). The space i{iu,v) : R ^ R x R},\\ ■ \\e) 
is clearly complete, and so {{4>{m), 4>(m),x)}m=o converges in this space. This yields 
= hm^^oo 0{m) (a;) e H'^''^{R) for all 1. Furthermore, we can take the 

limit m — > oo in ()19|) to obtain 

< Ci^[^](0), a;G[x„a;,], 

which shows that sup^-gj^.. ^.j ||0(a;)||i,2 is bounded. We obtain the pointwise bound by 
applying Theorem 4 to the sequence {0(m)}m=o taking the limit m — >• cxd. ■ 

Remark 5.1 The fact that the coefficients of (jH)) are independent of p means that the 
spatial derivatives 9^0 also satisfy this equation, and so Theorem 5 can be applied to 
these spatial derivatives. For example, if we assume that 0,ppp and (p^xpp are initially in 
(or more appropriately, that initially G H^'^ and (p^x ^ H"^'"^), then we can apply 
Theorem 5 separately to 0, and 0_pp to obtain solutions for which and its first 
two spatial derivatives are bounded on [a;j,a;c] x M. Likewise, if we specify data for (jHJ 
with smooth compact support, then the solution and all its spatial derivatives will be 
smooth, will have compact support on each slice x = constant and will be bounded on 
[xi,Xc] X R. 

4. Local energy measures 

In considering local observations of the energy content of the multipole field 0, two 
quantities are of relevance. These are the fiux of measured by an observer O moving 
on a timelike geodesic (with unit tangent u"): 

^ = m"9,0, 

and the "time-time" component of the energy-momentum tensor Tab of the field. For 
the massless scalar field, we have 

Tab = da4>db4> - ]^gabg'"^dc(t)dd(t), 



Bounds for scalar waves on self-similar naked- singularity backgrounds. 



13 



and the energy measured by O is 

Using the third part of Theorem 3 and an argument identical to that used to prove 
Theorem m we have global pointwise bounds on (xc — xY^'^cp^^- From this and the line 
element (jHl), it is straightforward to show that the term g'^'^dc4>dd4> is bounded on Ai. 
Thus finiteness of both S and JF follows if and only if the term u°'da4> is finite. In |3j , we 
showed that in the coordinates x, p, the components of the unit vector field tangent 
to O remain finite throughout A4 and (crucially) at Ti. Thus u"'da(f) is finite if and only 
if both of the terms (p^x and are finite. We know that the latter is finite, and so the 
local energy [S or JF) is finite at Ti. if and only if (p^x is finite at Ti. 

We can show that this term must in fact be finite by rewriting the wave equation 
(jHl) as a first order transport equation for (p^x with a source term involving only and 
its spatial derivatives. The fact that these terms are bounded enables us to write down 
a formal solution of this transport equation and hence demonstrate that (p^x is bounded 
in the limit as the Cauchy horizon is approached. So we define u := (p^x- Then (jH)) reads 

au^x + Wu,p + {a + 2(1 - n)(3)u = F{x, p), (20) 

where 

F{x, p) = -70,pp - {(3' + 2(1 - /€)7)0,p + [k{2 - k)^ + k(3' + ^{^ + l)e^]0. 
The characteristics of this first order equation are defined by 

^ = 2- 

dx a 

giving 

p=i+ r 2^ds=i+\{x), 

where ^ labels the different characteristics, and gives the value of p on a given 
characteristic as it intersects the initial surface x = Xi. We note that these characteristics 
are the outgoing radial null geodesies of the spacetime. The equation (pUj) can then be 
written as an ODE along individual characteristics: 

a-^ {u{x, ^ + A(x))} + {a + 2(1 - K)f3)u{x, ^ + \{x)) = F(x, ^ + A(x)). 



Defining 



J{x) = exp ^ I — i }R 



the solution of the ODE can be written as follows: 

J{x)u{x, ^ + X{x)) = u{x,, + r ^ + ^(^)) (21) 

Jx, a{s) 

We now consider the limiting value of u as the Cauchy horizon is approached, 
using the solution (PT|) . We consider first the case of smooth initial data with compact 
support. Noting that \imx-,xc A(a;) = +oo, we see that approaching x = Xc at constant p 
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entails ^ — > — oo. Thus for values of x sufficiently close to Xc, the characteristic through 
the point (x, p) meets the initial data surface x = Xj at an arbitrarily large, negative 
value of p, yielding u{xi,C,) = (where we have appealed to the fact that the data have 
compact support). 

Next, we note that it follows from Remark 15.11 and the boundedness of the 
coefficients of (j) and its spatial derivatives in F that F{x, p) is smooth, has compact 
support on each surface x =constant and is bounded on [xj, Xc] x M. The bounding term 
will be of the form Co£'[(^](0) + Ci-E[<^^p](0) + C2-E[<^^pp](0), where (as usual) the Ci are 
constants depending only on the functions G, ip and the angular mode number ^. Here, 
(p^p is to as is to (and likewise for the second derivative). Applying the mean 
value theorem for integrals gives 

r 44^(^' ^ + ^(^)) = ^ + ^(^*)) r 44 

Jx, a(s) Jx, a{s) 

for some x* G [xj, x]. 

From the definitions and (jlUj) of a and (3 and the assumptions of Section 2, we 

have 

a' + 2{l-K)(3 K + xlG\xc 



a 1 + x2G'(x, 



-1 + 0(1) 



Thus 



X — X^ + (J i , X — i> Xr 



K-1 



J(x) = Jo|x - Xcl + 0(|X - Xcl ), X^Xc 

where < Jo < oo and so 

J(x) a(s) ^' 



a; 



for some finite number Ji. 

Combining the results of the last three paragraphs yields the following: 

lim u{x,p) = F(x**,p)Ji (22) 
for some x** G [xi,Xc]. This proves the following result. 

Theorem 6 Let Xj G (0, Xc) anc? /et f,g & C^(M, M). T/ien t/ie unique solution 
(f) G C°°(M X [xj,Xc),M) o/ i/ie initial value problem consisting of (0) wt/i the initial 
conditions 4'\x=x, = f,4',x\x=xi = 9 satisfies limri:^xc4',x{x, p) G M for all p G M. 
Furthermore, this derivative satisfies an a priori bound of the form 

lim |0,.(x,p)| < CoEmO) + CiE[0J{O) + C2E[0^J{O), p G M. (23) 



Again, invoking the density of C^(]R) in certain Sobolev spaces, we can obtain a 
more interesting result where the field (p does not vanish at the origin. 
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Theorem 7 Let Xi E (0,Xc) and let f G H^''^(M.), g E if^'^(]R). Then there is a unique 
solution (f) E C([a;j, Xc], -ff^'^(M)) of the initial value problem consisting of (0j with the 
initial conditions (f)\x=xi = f,4',x\x=xi = 9- This solution satisfies 

lim |0,.(x,p)| < CoE[if]{0) + C,E[ipJiO) + C2^[^,J(0), p G R. 

Proof: The existence and uniqueness part of the resuh is an apphcation of Theorem El 
and Remark l5.ll We require the third derivative of / to be in to ensure finiteness of 
E[ip^pp] (0). The a priori bound obtains by applying TheoremlHlto a sequence of solutions 
{4>{m)}m=o of dHl) which Satisfy 4>(m){xi, p) ^m^oo f in if^'^(M) and (f>(m),x{xi, p) ^m^oo 9 
in i7^'^(]R). The bound (j2SI) applies to each member of this sequence, and so applies to 
the solution in the limit. ■ 



5. Comment on massive fields and non-minimal coupling 

We consider briefly a massive minimally coupled fleld 0, satisfying the Klein-Gordon 
equation 00 — = 0. The mass parameter m > introduces a length scale, and so 
must break self-similarity. This is reflected in the presence of a radially dependent term 
in the wave equation: the massive fleld satisfles an equation differing from ((7j) only in 
the term in 0: 

a0,xx + {m^e'^P + i{i+l))e^(j) = 0. (24) 

The presence of this additional term invalidates some but not all of the results above. 
We are able to retain the results relating to wave packets, i.e. results relating to the 
case where the initial data are smooth with compact support. In this case. Theorem |21 
applies. An amended version of Corollary 01 also applies: we can obtain a bound of the 
form 

where p((p(0)) is the Lebesgue measure of the support of the initial data <p(0). This 
number will determine - via the characteristics of the system (fT^ - a maximum for p on 
the support of the solution at time t, and so will determine a bound for |-B(t)| that 
will depend on t and p((p(0)). This will have the form indicated above. Noting that 
Lemma El remains valid for m 7^ 0, we can conclude that Theorems El and 01 hold for a 
massive, minimally coupled scalar fleld. However the arguments used to prove Theorem 
5 break down due to the presence of the Lebesgue measure of the support of the initial 
data in the bound above. It would be of interest to determine if another approach could 
be used to obtain the equivalent of these results in the massive case. This should be 
feasible, especially as what one is most interested in doing is extending the function 
space at the origin p = —00, where the radially dependent term is exponentially small. 
The argument of Theorem 6 also breaks down, as we will not have a global bound for 
the function corresponding to F. 
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We could also consider other couplings, so that satisfies e.g. 00 — = 0, where 
^ is a constant and R is the Ricci scalar. Here, an additional term of the form q{x)(j) 
is present in the wave equation, so self-similarity is preserved (as expected). This will 
make a crucial difference in Lemma and it does not seem possible to draw any general 
conclusions as to the continued applicability of this lemma, and hence of the validity 
any of our principal results. We expect that similar results could however be obtained 
on a case by case basis, where the background geometry is specified. 

6. Conclusions 

We have shown that in self-similar collapse to a naked singularity, the multipoles of 
a massless scalar field propagating on the background spacetime remain finite as it 
impinges on the Cauchy horizon. This has been shown to be true of different measures 
of the multipole field: its point- values, its norm and different local {S,J-') and global 
(E) measures of its energy content. We have concentrated on the field = ^e"'', and 
so it is instructive to consider the implications for of (in particular) the pointwise 
bounds obtained. We have shown that r'^(j){x,r) and r'^0 a.(a;, r) are bounded on the 
Cauchy horizon for all k G [0, k*) where > is determined by the background 
geometry. Thus the values of and (p^^ on the Cauchy horizon can only diverge at the 
singularity. In particular, any possible such blow up does not make itself felt along the 
Cauchy horizon: it is confined to the central singularity. Similar (but more limited) 
results apply also for massive scalar fields. The question of whether or not these results 
will also apply to the full field obtained by resumming the multipoles remains open, 
but it should be noted that the multipoles themselves have an independent physical 
significance. For example, one would expect that the scalar radiation field is dominated 
by the i = contribution. 

We interpret our results as providing evidence that these naked singularities may 
be stable under perturbations of the background spacetime, and so may constitute a 
challenge to the cosmic censorship conjecture. In particular, we consider that these 
results add weight to the existing evidence that GRLP spacetimes admit stable naked 
singularities 0. As there are no unstable modes of perturbations impinging on the 
singularity, it is likely that these spacetimes give rise to the type of initial data studied 
here. This is also possible for the single-unstable mode critical spacetimes ^T]. In these 
spacetimes, the single unstable mode shows a characteristic divergence as r ^ along 
surfaces x = constant. If this divergence is sufficiently mild, this mode would correspond 
to initial data of the type studied here, where weighting by a power of r brings the data 
into a certain Sobolev space. The situation described in the previous paragraph could 
then hold, where the divergence at the singularity is not felt along the Cauchy horizon, 
yielding stability of the horizon. 

Let us suppose then that the results found here for a minimally coupled scalar 
field carry over to the metric and matter perturbations of the background spacetime. 
This would indicate only linear stability of the naked singularity. An entirely different 
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picture may emerge when full (non-linear) stability is considered. For example, in the 
case where the background energy momentum tensor is that of a minimally coupled 
massless scalar field (Einstein-Klein-Gordon spacetime), examples of naked singularities 
are present in the self-similar case jT21[ISI- Our results apply to these examples, and so 
provide evidence of linear stability. However these have been proven to be non-linear 
unstable ^3]. Thus our analysis is probably best interpreted as demonstrating the 
validity of a necessary - but nonetheless nontrivial - condition for stability. 

Finally, we note the connection with generalised hyperbolicity |15j /wave regularity 
[ini HZj /quantum regularity [IHl dl 1201 of naked singularities. These are broadly 
similar ways of characterising whether or not a naked singularity does in fact provide 
a fundamental barrier to predictability in spacetime. As a next step in the study 
of naked singularities in self-similar spacetimes, it would be of interest to determine 
whether or not these singularities pass the tests of generalised hyperbolicity/wave 
regularity/quantum regularity. A minimum requirement for passing such tests is that 
the field remains intact as it impinges on the Cauchy horizon. As shown above, this 
requirement holds. The more complicated question of the well-posedness of the field to 
the future of the Cauchy horizon remains to be answered. 
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